Pseudospin entanglement and Bell test in graphene 
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We propose a way of producing and detecting pseudospin entanglement between electrons and 
holes in graphene. Electron-hole pairs are produced by a fluctuating potential and their entangle- 
ment is demonstrated by a current correlation measurement. The chirality of electrons in graphene 
facilitates a well-controlled Bell test with (pseudo-)spin projection angles defined in real space. 
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The entanglement between internal degrees of freedom 
of an electron and a hole in the Fermi sea is of both, fun- 
damental and practical interest. It has been recognized 
as a form of entanglement that does not require many- 
body interactions [l| and comparatively simple ways of 
generating it experimentally have been proposed [l|, Q ■ 
Yet, the experimental demonstration of electron-hole en- 
tanglement in solid-state structures is still outstanding. 

Very recently Neder et al. [3] have implemented a 
quantum Hall interferometer that had been proposed |2| 
for the generation and detection of electron-hole entan- 
glement. This interferometer has allowed the observation 
of the "two-particle interference" that is at the core of 
electron- hole entanglement 0]. While this is a first indi- 
cation of entanglement production in the experiment of 
Neder et at, a conclusive verification of that entangle- 
ment has not yet been achieved. It has been proposed 
P, 0, 0, m that the generated entanglement is best ver- 
ified through the violation of a Bell inequality. Experi- 
mentally, the demonstration of such a violation, however, 
meets with significant challenges. First, the decoherence 
in the interferometer of Ref. [3] needs to be reduced signif- 
icantly in order to safely preserve the entanglement from 
the time of its production to its detection. In addition, 
a test of Bell inequalities requires measurements of spin- 
1/2 degrees of freedom along variable quantization axes. 
In the setup of Ref. ^ these quantization axes are defined 
by scattering amplitudes that are poorly controlled ex- 
perimentally, requiring an implementation through trial 
and error. 

In this Letter we develop a way of generating and de- 
tecting electron-hole entanglement that does not suffer 
from the above mentioned problems. We propose to en- 
tangle the pseudospin f?\ degree of freedom of electrons 
and holes in graphene 8, 9, 10. by means of the "pump- 
ing" mechanism of Refs. [ul, liij. The typical energy 
scales in graphene are considerably higher than those in 
GaAs, which has for instance allowed an observation of 
the quantum Hall effect at room temperature [l3|- At 
comparable temperatures one thus expects the decoher- 
ence in graphene to be much weaker than in the experi- 
ment of Ref. addressing the first of the above issues. 
In addition, we formulate a Bell test through current 



correlation measurements that overcomes the mentioned 
problems ofpreviously pursued entanglement detection 
schemes [H, Electrons in graphene have a definite 

chirality (for a certain bandstructure "valley"), moving 
in the direction of their pseudospin. The pseudospin of 
an excitation can thus be measured through its direction 
of motion. This affords a Bell test with straightforward 
and transparent control of the (pseudo)spin quantization 
axes that are now defined in real space. The Bell test 
proposed in Refs. [l|,S,lli] is only valid in the regime of 
temperatures T that are low compared to the voltage V 
applied to the interferometer: kT ^ eV. Its application 
at finite temperatures faces a problem that has been dis- 
cussed recently in Ref. fij . The authors of Ref. [ij suggest 
a cure of that issue whose experimental implementation, 
however, is challenging: It requires the addition of res- 
onant levels to the setup. Here we avoid the problem 
pointed out in Ref. ^] by a suitable postselection of the 
entangled electron-hole pairs. That selection is imple- 
mented simply by subtracting the thermal background 
from all measured current correlators. 

Setup: We consider a sheet of ballistic graphene at low 
temperature T and nonzero Fermi energy ep ^ kT in a 
vanishing magnetic field B. The sheet is well-coupled to 
an electron reservoir along its rim, as shown in Fig. [T] 
We formulate the low-energy Hamiltonian of graphene in 
single- valley form through a unitary transformation that 
renders the Dirac model valley-isotropic 



Ho 



V cr ■ p. 



(1) 



Here, cr = (axyTy) is a vector of Pauli matrices in pseu- 
dospin space, p is the electron momentum, and v is the 
Fermi velocity. The graphene sheet is subject to a local 
fluctuating potential, as described by the Hamiltonian 



i/ex(i) = / dxu{x) eFexW i^Hx) ■ i^{x). 



(2) 



Here we have switched to second-quantized notation with 
electron annihilation operators tp that are vectors in pseu- 
dospin space. The fluctuating potential is focused on a 
small region of spatial extent l^x in the middle of the 
graphene sheet and it is centered around the origin of 
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FIG. 1: Proposed setup: A localized fluctuating electric po- 
tential Vex produces entangled electron-hole pairs in the cen- 
ter of a graphene sheet. The excitations are either drained 
at the rim of the sheet or leave into tunnel contacts at lo- 
cations Tea (where a takes the values ±a, ±a', ±6, or ±&'). 
The generated entanglement has signatures in correlations be- 
tween the currents into these tunnel contacts. The starlike 
setup that is shown allows the maximal violation of a Bell 
inequality. 



our coordinate system. Its shape is given by a func- 
tion u that is normalized to fcp Jdxu{x) = 1 [for in- 
stance u{x) = exp(— |a;p/2Z^x)/27rfcpZ^^]. The length Zox 
is assumed to be large compared to the lattice spac- 
ing ^lattice, but small compared to the Fermi wavelength 
27r/fcF, that is Ziattico ^ 'ox ^ 27r/fcF. We assume 
that the frequency spectrum of the potential correlator 
cy(— — J dt eyip{—iujt)(Vcx{t)Vc^{0)) is relatively flat 
until it vanishes above a high-frequency cut-off w « 
with fcT < 17 < £f [UHllH. The fluctuating poten- 
tial Vex creates pairs of electrons and holes (in the sense 
of an electron missing in an otherwise filled Fermi sea) 
that propagate outward before they are reflectionlessly 
drained by the reservoir surrounding the graphene sheet. 
On their way to the rim they are able to leave through 
tunnel contacts a with sizes la ^ into additional elec- 
tron reservoirs at locations where |q;| = 1 and 
Ta » 27r//cF (see Fig.[T|), 



dx ^\x) ■ 



w^{x) 
w^{x) 



^7 



h.c. 



(3) 



Here, the functions Wa are centered around 
and the operator annihilates electrons in the reser- 
voir of contact a. All electron reservoirs are in thermal 



equilibrium with the graphene sheet. Every tunnel con- 
tact a has one counterpart —a in direction —a. 

Entanglement production: We first consider the ex- 
citations created by a short potential pulse, eVc-x{t) = 
C,5{t — tcx). In first-quantized form the low energy con- 
tribution {\p — p'\ <C fcp) to the electron-hole pair that is 
produced at first order in C reads (we set ?i = 1) 



\m)c 



pp' 



ur'ii)'^), 

(4) 

where "f and J, specify the pseudospin direction and |)°' 
and 1)'^ are electron and hole amplitudes, respectively 
(with the convention that a hole has the same pseudospin 
as the electron that it replaces). The pseudospins of the 
electron and the hole described by |V')ci Eq. (jl]), are en- 
tangled. They form a so-called Bell pair. A source with 
a periodically varying potential Vox(t) serves as a steady 
supply of such Bell pairs. Excitations that appear at 
higher order in ^ are negligible if |eI4x(OI ^ for all t, 
which we assume henceforth. 

Entanglement detection: The Heisenberg equations of 
motion corresponding to the Hamiltonian Hq, Eq. ([1]), 



p — 0, CT — 2vp X (T, X ~ va, 



(5) 



show that the pseudospin of an electron is not con- 
served. We therefore postselect orbital states of the form 
\p)a — {\p(^) + I ^pct))/V^. The initial wavefunction 
\ip)c, Eq. (j4|), factorizes into an isotropic orbital part and 
a pseudospin part. It follows that at t = fox also the 
pseudospin state of all postselected electron-hole pairs is 
given by the pseudospin part of Eq. ([4]) and entangled. 
Moreover, the (unnormalized) density matrix of a post- 
selected electron p^^, = {p\aP\p')a: P being the density 
matrix of the electron before postselection, takes the form 
p^p,{t) = cos[2wp(t-iex)]cos[2V(t-iox)]Ppp-(tex). After 
normalization is time-independent. The pseudospins 
of the postselected excitations are thus conserved and so 
is their entanglement. We propose to verify that entan- 
glement by violation of a Bell inequality. This requires a 
measurement of the postselected pseudospins with vari- 
able quantization axes. The tunnel contacts a serve that 
purpose. To see this we integrate Eqs. (O to find 



X{t) = X{tcyc)+v{t-tcx) 



[P(<cx) 



cr(icx)]p(icx) , / 1 



\Pito.W 



(6) 

Consider an electron (before postselection) that is pro- 
duced aX t = iex and (a;(tox)) = 0, such that semiclassi- 
cally the first term in Eq. ^ vanishes. The third, oscil- 
latory term in Eq. ([6]) is smaller than the second one by a 
factor (fepTa)"^ and negligible in our limit. Assume that 
this electron is detected in contact a at time t, such that 
x{t) = raOi. This projects the initial state of the electron 
onto eigenstates of x{t), Eq. ([S]), with eigenvalue raOt. 
To our accuracy, when only the second term in Eq. ([5]) is 
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FIG. 2: Pseudospin measurement: Electron and hole excita- 
tions are generated by \4x at a; = 0. They contribute to the 
current into contact a at a; = r^a if their velocity sat- 
isfies ilia. The pseudospin cr of these excitations is fixed 
due to their chirality: one has <t 1 1 i and cr thus points "up" 
along the quantization axis a. The opposing contact —a at 
X = —raO. collects excitations with pseudospin "down" along 
the same axis. 



relevant, such a measurement projects onto amplitudes 
withpllo:. Moreover, substituting p(tox)/|p(^cx)| = 
into Eq. ([6]) (the sign of p{tc^) is arbitrary since in enters 
quadratically), we see that a measurement of an electron 
in contact a projects onto amplitudes with cr{tc^) -o; 1. 
Likewise, detection of an electron in contact —a projects 
onto amplitudes with (t(<ox) • a = — 1 We conclude 
that the tunnel contacts a collect currents la of electrons 
and holes with a definite initial pseudospin when mea- 
sured along quantization axes a that are defined by the 
locations of the contacts a in real space, as illustrated in 
Fig. [2] As discussed above, the initial pseudospin state of 
the created electron-hole pairs before postselection equals 
the pseudospin state of any of the postselected electron- 
hole pairs. Also the pseudospin of the postselected pseu- 
dospins may thus be inferred from a measurement of the 
currents /„. 

In order to demonstrate the entanglement of the post- 
selected electron-hole pairs we formulate a slightly mod- 
ified Clauser-Horne-Shimony-Holt inequality 



B < 2 with B = Cab + Ca'b + Cab' " Ca', 

in terms of symmetrized correlators 



(7) 



Cfc = i (^1 (d . <T°') (5 . ^ . ^ci) (^a-a^m (8) 

of an electron pseudospin cr''' and a hole pseudospin cr^, 
projected onto the unit vectors a and b. Only entangled 
electron-hole pairs can violate the inequality ^ and the 
parameter B is thus an "entanglement witness" [l^]. 

At zero temperature and in our limit of dilute electron- 
hole pairs |eI4x| ^ ef one shows along the hues of Refs. 
0, Q that the parameter Cab can be expressed through 
zero- frequency current correlators. One needs to cor- 
relate the currents of electrons and holes with definite 
pseudospins when measured along the quantization axes 
a and b, respectively. The reasoning of the paragraph 
around Eq. ^ allows us to relate these currents to the 



tunnel currents I±a and I±b into the contacts a = ±a 
and a = ±b. In case all tunnel contacts couple with 
the same strength to the relevant electrons or holes one 
concludes in this way that 0, ^ 



Cab 



(9) 



where Ca.b = j dt (SIa{t)SIb{0)) has to be substituted for 
Ca,b- Unlike in the proposals of Refs. here not every 

produced electron-hole pair is detected. To find Ca.b hi 
the general case, when the tunnel coupling strengths to 
different reservoirs a are not equal, one therefore has to 
normalize the correlators Ca,b by the detection probabil- 
ities Wa, 



Ca,b = {WaWby J dt [{Sla{t)dlb) - {SIa{t)dIb)v^^=o] ■ 

(10) 

The probabilities Wa can be measured through the AC- 
response of the currents to the excitation potential Vox 
at frequencies a; <C min{t)//Q,, ep}, 



Wa = 



e2Kx(w) 



(11) 



Extra care has to be taken at finite temperature. 
We avoid the issue pointed out in Ref. Q here by an 
additional postselection of the electron-hole pairs for 
which Eq. ([7]) is evaluated (see Appendix A for the de- 
tails). This selection is made by the subtraction of the 
equilibrium current correlations in Eq. (|10p . We show 
in Appendix A that after that subtraction a violation 
of the Bell inequality ([7]) is an unambiguous signature 
of entanglement also at finite temperature, as long as 
v/\raa±rtb\ < fcT < fi. 

Predictions: In our limit of small Wa (that is tunneling 
contacts), v/\raa ± rbb\ <C A:T <C <C £f, leKxl ^ ep, 
iex ^ 27r/A:F, and la ^ ra we find (see Appendix B) 



Ca,h = ^3-^(1 + d • b) / de de'cvie-e'). (12) 

We conclude that the pseudospin correlators measured 
through Eqs. ([9]) and pO|) take the form 



Cab = a - b, 



(13) 



which is immediately shown to violate the modified 
CHSH inequality ([7]) for appropriate choices of the vec- 
tors a, a , b, and b . The maximal violation B = 2^/2 for 
instance can be achieved in the symmetric starlike setup 
shown in Fig. [1] where the vectors d , 6, d, and b are 
separated by successive 45° angles. 

For Eq. ([13]) to hold the created excitations must not 
change their pseudospin on the way to the tunnel con- 
tacts, for instance through decoherence. We thus as- 
sume excitation energies such that the inelastic mean free 
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path Zin is long, ^ ra- This condition is fulfilled at 
n <C mm{vpiikF , £f /VkvrZ} , where Wph is the phonon 
velocity in graphene 1^, 2^. As anticipated, for suitable 
parameter values decoherence through inelastic processes 
in our proposal is already suppressed at temperatures 
much higher than those in the experiment of Ref. 
We need in addition that also the elastic mean free path 
lei is long, Zei ^ fa- Long elastic mean free paths have 
been found in suspended sheets of graphene [2l| . We con- 
clude that the setup depicted in Fig. [Hallows to generate 
and conclusively demonstrate electron-hole entanglement 
if min{z;phA:F , Ef / V^pf q } > il > fcT > v/\raa ± rbb\. 

Discussion: We have proposed a way of creating and 
verifying pseudospin entanglement of electron-hole pairs 
in graphene. Bell pairs are produced by a fluctuating 
potential and their entanglement is demonstrated after 
postselection through violation of a Bell inequality. The 
quantization axes in the requisite pseudospin measure- 
ment are defined by the locations of tunnel contacts in 
real space. This simplicity of the pseudospin measure- 
ment is bought at a price: the postselected Bell pairs are 
not easily separated spatially, since the advocated pseu- 
dospin measurement is nonlocal, as shown in Fig. [21 The 
produced Bell pairs, however, are entangled also with re- 
spect to their intrinsic spins l23|jwhich entanglement is 
readily spatially separated [llj, [l^]. The proposed ex- 
periment is thus an intermediate step toward the gener- 
ation and manipulation of spatially separated Bell pairs 
in electronic structures. Entanglement is generated by 
a mechanism that is able to produce spatially separated 
Bell pairs, but it is detected before that spatial separa- 
tion is achieved . The proposed detection mechanism 
affords three major advantages: i) it suffers less from 
decoherence than previously pursued implementations of 
particle- hole entanglement; ii) it allows a well-controlled 
Bell test with clearly defined (pseudo)spin quantization 
axes; iii) it avoids problems with earlier proposals of 
electron-hole entanglement detection at finite tempera- 
ture. Our proposal thus overcomes some critical hurdles 
on the way to an observation of particlc-hole entangle- 
ment in electronic structures. 

The author thanks W. A. de Heer, P. N. First, and L. 
You very much for discussions. 
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APPENDIX A: BELL TEST 

The relation between current correlations and the 
(pseudo)spin correlator Cab, Eq. ([9|), has been proven 
m Refs. BM 

at zero temperature in the limit of dilute 
electron-hole pairs (in the system discussed here that is 
|eX4x| ^ £f)- It has been pointed out recently [q] that 
extra care has to be taken at finite temperature. Ther- 
mal excitations in the reservoirs that collect the mea- 
sured currents are then able to generate electron and 
hole currents that flow in the "wrong" direction: from 
the reservoirs into the system that contains the entan- 
gled electron-hole pairs. These currents pose a serious 
problem and they can invalidate a Bell test along the 
lines of Rcfs. [4, 5] at finite temperature. In Ref. an 
alternative, energy-selective detection scheme has been 
proposed that does not suffer from the same problem. 
Experimental implementations of energy-resolved detec- 
tion, for instance through resonant levels [2^, are, how- 
ever, expensive. In contrast to the situation studied in 
Ref. the entangled electrons and holes in our proposal 
differ in energy. This affords a simpler cure of the prob- 
lem. 

The typical energy separation between the electron and 
the hole in the produced Bell pairs is fi. This allows 
us to restrict our attention by means of postselection to 



electron-hole pairs with electron energies e > ep + w and 
hole energies e < ep — u), where we choose kT <^ uj fl. 
Thermally activated electron-hole pairs in the graphene 
sheet and the reservoirs coupled to it have typical en- 
ergies |£ — EfI, |e — £f| ~ kT. All but an exponentially 
suppressed number of them are excluded by the above 
postselection. This avoids the problem that has been 
pointed out in Ref. At the same time our postselec- 
tion includes almost all excitations created by Kx, that 
have typical energies |e — ep], |£ — SfI ~ ^- The above 
postselection may thus be implemented approximately by 
subtracting the statistical (thermal) contributions from 
all measured current correlators, leaving only correlations 
due to Veic- These thermal correlations, in turn, may be 
inferred from a measurement of the respective correla- 
tors in equilibrium, at Vex = 0. The above postselection 
should thus solve the problem pointed out in Ref. @ 
with very moderate additional experimental effort: it is 
implemented by a second measurement in thermal equi- 
librium, as expressed in Eq. (fTO|) . Below we prove this 
expectation correct. 

The two-particle density matrix after our postselec- 
tion of electron-hole pairs in the conduction band with 
symmetrized orbital states and \p)^ for electron and 
hole, respectively, at electron momenta p > (ep +lu)/v 
and hole momenta p < (ep — lj)/v reads 



pf^^,%?[p,p]p ,p) oc 6(wp - ep - w)9(ep - w - i;p)e(wp' - ep - w)6(£p - w - up') (14) 

X {p\i {af {p\i (^iVcondbx' wT \p')i 

\Q{x) = 1 for X > Q and 8(2:) = otherwise], where the pseudospins of the electron and the hole are denoted a and 
(T, respectively. Pcond is the density matrix of all electron-hole pairs in the conduction band of the graphene sheet 
(before postselection) . Of all momentum-restricted electron-hole pairs we then select those whose electron and hole 
are found within a radial distance |r — f| ^ Ar of each other (in real space), as described by the reduced density 
matrix 

« rdrdfe-(^-'^'/'^^'pfr^^^^{r,f;r,f). (15) 
Jo 



Here, p**'"^ is the Fourier transform of p®^"'' [with the 
sign convention of Eq. pS]) ] and R is the radius of the 
graphene sheet. As explained in the main text, p°^-^°'°' 
is independent of a. and a for the state Eq. We 
have therefore suppressed these indices of p^^™. For the 
proposed Bell test one needs to measure the pseudospin 
correlators Cab, Eq. ([8]), evaluated for the postselected 
electron- hole pairs. They follow from 

cl = J Trp^P'"(l + d . ^-•)(1 + 6 . ^1^) (16) 

as Cab = cl + cl - - Cl_^ - Cl_fj - 



C-6a + C'-a-fc + C'-fc-a- Wc assumc > V / Ar > 
max{n{eVex/eF)^,kT{eF/eVox)^{kT/n) exp{-2cu / kT)} . 
In this limit it is very unlikely for two holes to be within 
a distance Ar of the same electron (or vice versa). The 
correlators Cl may then be measured by "coincidence 
detection" [J 0], correlating momentum-projected 
electron and hole pseudospin densities, 

cl (X r drdfe-^^-'^'/^'^'^^ {nl^\r)nl''{f)). (17) 

Here, the densities n^^"' and n)^ are defined in terms of 
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the momentum-projected electron annihilation operators The correlator Cjj^, however, may be accessed also 



£f) - w] ^ Vjjai'Tjpa 
V=±l 

(18) 

as n'j/{r) = [ip^^{r) ■ vl^][vl} ■ '0q(7')]- The index /i takes 
the values el or h and we have introduced 7°' — 1 and 
7*^ = —1. The vectors v\ are the spinors corresponding to 
pseudospin up along the directions a and the pseudospin 
matrix Vp projects onto the conduction band. Note that 
at this point the choice of momentum direction of the 
postselected excitations is arbitrary: in the state ^ am- 
plitudes with any momentum have the same pseudospin. 
Our above choice of momenta along the desired pseu- 
dospin quantization axis a has been made merely for 
convenience. In our limit, when all electron- hole pairs 
are well-separated from each other (much farther than 
Ar), the density correlator of Eq. p7|) may be replaced 
by its irreducible contribution 



cx / drdr e 




-(r-r)^/2Ar^ 



{5nl%)5nl'^{f)) (19) 



0, The momentum-projected density is in princi- 
ple experimentally accessible, for instance by tunneling 
electrons and holes into additional reservoirs that couple 
symmetrically to two contacts a and —a via momentum- 
dependent tunnel amplitudes and w^. Pseudospin- 
resolved amplitudes , as described by the tunneling 



Hamiltonian 



iJ^^ with 



h.c. 



then allow to obtain the correlator C|{, as 
from current correlators 



(20) 



(21) 



(22) 



where 



(23) 

We have introduced one reservoir for every pseudospin 
state, with corresponding electron annihilation operators 
^res/i^ The amplitudes are peaked at the radii r^. 
In the step from Eq. (fT9|) to ((2T|) we have assumed time- 
translational invariance and r°' = . Note that the error 
due to thermal excitations flowing from the reservoirs 
into the graphene sheet (the origin of the problem pointed 
out in Ref. |6i]) is here exponentially suppressed. It is of 
order exp(— w/fcT) since it is only pairs of excitations 
that differ in energy by at least uj that contribute to the 
correlator c^®'"*^. 



without energy-selectivity, 
/o°° if exp(ipr) E„=±i 'Pv^^. 



through measurements of correlators 

cL = j dte-^^'^''^^^\6ll{t)5llm (24) 
of currents 

Il^iel drw^{r) [^Ur] ■ ^^K^t . ^-) + h.c. (25) 

where ^a{f) = 
ripoL (and a correspond- 
ing tunneling Hamiltonian Hj.^). This is seen easiest 
by separating two contributions to the above current 
correlators from each other: First there are contribu- 
tions due to the electron- hole pairs created by T^x- 
We denote these contributions to c^°'~^ and by a 
subscript "ex", cj^'~^ and cj^, respectively. Second, 
there are contributions from statistical correlations 
of excitations due to the Pauli principle, denoted 
by a subscript "stat," cJ^^tT'^ ^-i^d (^lta,f We have 
= + cLt and likewise c^''^-^ = cj^'-i^ + clf~^. 
The typical energy separation between the electrons 
and the holes created by Vox is vp — vp « Q. For 
the currents /Jxa carried by these excitations we thus 
may approximate /J^^ = /J^'q, + /J^a + 0{uj/Q) and 



J 



■ Cg°'j,J^ + 0{uj/fl). Under the assumption 
kT 3> v/\r^azLr^b\ made in the main text the statistical 
correlations c^^^Lt identical to the equilibrium corre- 
lations, cLt = cT|k.=o + 0[/cTexp(-A:T|<a ± r^b\/v)] 
(see Appendix B). Moreover, statistical fluctuations do 
not contribute to the momentum projected correlator 
by our definition of V-^: OV'avC^'- = 0. 

We conclude that to leading order in our limit 
the momentum projected irreducible current cor- 
relator c^°'~'^ may be expressed through the cor- 
responding correlator without momentum projec- 
tion after subtraction of its statistical background, 

Tol-h Tcl-h _ Tol-h , Tcl-h _ T _ T T 



such that 



ba 



c 



"ab I Vo, 



Vcx=0 



(26) 



The pseudospin-resolved tunneling amplitudes assumed 
above are rather unrealistic. Alternatively the pseu- 
dospin currents can be measured by introducing sepa- 
rate reservoirs for the contacts a and —a, as explained af- 
ter Eq. ([6]) of the main text. One further shows straight- 
forwardly along the lines of Refs. [1,0] that the correla- 
tors c, Eq. ((24)) may be replaced by zero- frequency cor- 
relators [with an error of 0{v/nAr)]. Eqs. dH) with ^ 
then follow after a normalization of along the lines of 
Refs. [3, 0|. Values for uj and Ar that satisfy all of the 
above conditions can be found provided that kT <C A 
violation of the inequality ([7]) with the correlators Eqs. 



7 



([9]) and ((TO)) is thus indeed an entanglement witness to 
leading order in ?;/|r^a ± r|^b| <s:. kT fl. 

Strictly speaking, a Bell test has to demonstrate that a 
violation of Eq. ([7]) found through the current correlation 
measurement described above is due to the entanglement 
of p'^P'" rather than the deviations of the true pseudospin 



correlators from Eq. ([9]) that appear at nonzero tempera- 
ture (even though those are suppressed in our limit). Col- 
lecting all the sources of such deviations mentioned above 
we find that the Bell parameter B measured through Eq. 
© is related to the Bell parameter ^B'^p'" corresponding 
to the actual pseudospin correlators Eqs. ^ and as 



il ilAr V \ sf 



34 ■ 



kTAr kT f ep 



-2w/fcT 



EkT ( ep 



(27) 



up to terms of higher order in small quantities, with positive constants gj and gaa' that are of order unity. The 
summation in Eq. (|27p runs over all measured combinations of tunnel contacts a, a'. Entanglement is conclusively 
demonstrated if one finds a violation of the inequality S'^p'" < 2, that is if 



I — 



A direct violation of Eq. ((28|) requires knowledge of the 
parameters gj and gaa'- They can in principle be deter- 
mined if the frequency spectrum of the voltage correlator 
cv is known. More practically, however, one may mea- 



sure the dependence of B on parameters such as T and 
Kx and establish an inconsistency with Eq. pS)) for a 
particular choice of uj and Ar. For example, one may 
choose (jj = ^kTVL and Ar — e^v /VLeVo^^ such that 



B = B^^'-{Vkfn,eFv/neV^.)+gi^+{g2+93)'-^+9, (l^) ^ (^) ^ ^"'^^ + 



I fcT / £F \ -kT\ra&-r^,&'\/v 

aa' ^ ' 

An analysis of p^P'", Eqs. with (|15p . shows for the same choice of uj and Ar 



S''P'"(\/fc™,eFf/f^eV;x) = S"P'"(vfc7^,eFw/^2el^o) +.91 



-53 



■54 



^2 / VeV^o 



(29) 



(30) 



up to terms of higher order in small quantities, with 
positive constants g^ that are of order unity. Suppose 
that with this choice of lo and Ar one measures a pa- 
rameter S — 2 > that varies only by a small frac- 
tion ABjiB — 2) <C 1 over a range of temperatures 
T e [ro,2To] and potentials l^oxl € [Vb,2Vb] (Vq > 0). 
Eqs. ((29|) with (j30p prove this measurement result to be 
inconsistent with the assumption that no entanglement 



is present in the electron-hole pairs measured at Tq and 
Fq. It excludes that S"P'"(yHbn, eFf/f^eFo) < 2 (note 
that according to Eqs. and all T- and Vox- 

independent contributions to B have a negative sign). 
The result S - 2 > with variation AS/(S- 2) < 1 over 
the above parameter range, however, is predicted to be 
found in the presence of entanglement (for appropriate 
choices of the pseudospin projection directions such as 
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shown in Fig. [T]). In such a measurement the generated H^ a- Lowest order perturbation theory in iS and 14, 

entanglement can thus be verified conclusively. results in 



APPENDIX B: CURRENT CORRELATORS 



In this Appendix we obtain the correlators Cab 
J dt {SIait)6IbiO)) from the Hamiltonian H = Hq + Hc: 



^afe 



e"' /■ de de' , _ , , 

-— / dXadx'^dxbdx'f,——TTg{x'i,-Xa,s)[f'',Gais)]Wa{xa,x'Jg{x'^~Xb,e')[f'',Gbie)]Wbi^^ (31) 

Z J ZTT Zn 



describing statistical correlations. Here, g and Ga, are the Green functions of electrons in the graphene sheet and 
in reservoir a, respectively. They are matrices with Keldysh- 27| and pseudospin indices. We expand the Green 
functions g{x) asymptotically assuming k-p\x\ ^ 1 [28| . as it is appropriate in our limit kpra ^ 1- The matrix is 
the third Pauli matrix in Keldysh space and the unit matrix in pseudospin space, while Wa is unity in Keldysh space 
and 



Wa (Xa , X' 



<(a;aX*«) w^ixa)w^*{x'j 
wl{xa)wt*{x'a) w^{Xa)w^*{x'J 



(32) 



in pseudospin space. 

At the next to leading (second) order in Vox we find a contribution 



f de de' 

-— / dxadx[^dxbdx'i,dxc^dx'^^——cv{e-e')uixc^)u{x'^j 
Z J Zn Zn 



(33) 



X Tigix'b - a;cx,e)T^g(a;cx - Xa,e')[f'',Ga{e')]Wa{xa,x'^)g{x'^ - x'^^, e')f'' g{x'^^ ~ Xb,e)[f'',Gbie)]Wb{xb,Xb) 
due to electron-hole pairs excited by the fluctuating potential Vex and a contribution 

^^2)stat ^ J ^^t^ ^^/^ ^^^^ ^^/^^ cv(e - e')u{xcx)u{x'^^) (34) 

X Tr {g{x'b - Xcy^e)^ g{xcy, - x'^^,e')f^ g{x'^^ - Xa,e)[f^ ,Ga{e)]Wa{xa,x'^)g{x'^ - Xb,e)[f'' ,Gb{e)]Wb{xb,x'b) 
+ 9{x'b - Xa,e)[f'- ,Ga(.£)\Wa{xa,x'^)g{x'^ - x^y,, e)f'' g{xex - x'^^,e')f^g{x'^^ - Xb,e)[f'' ,Gb{e)]Wb(xb,x'b)] 



which describes statistical correlations due to a renor- 
malization of the tunneling amplitudes w by V^x- 

To lowest order in our limit <C £f and /ox, <C \ra\ 
we have w'^[x) = {v / \/kF)Wa5{x — rao) and u{x) = 
S{x)/kp. Eq. ([33]) then evaluates to 



c'iT = ^i^ + a-b)WaWb I de 



with the tunneling probabilities 



de' cv{e — e' 



(35) 



r 



where = * In ("w^a '^q*/ I'w^a |) and Va is the Fermi 
velocity in reservoir a. The above perturbation expan- 
sion in w is justified if ?D <C 1. 

The contribution c'^'*^**^', Eq. is exponentially 

suppressed as exp(— fcT|raa±7'h6|/i;) by thermal dephas- 
ing in the limit kT ^ v/\raa ± rbb\ taken in the main 
text. After the subtraction of the statistical fluctuations 
c^^^ and the normalization performed in Eq. (|10[) the cor- 
relator Ca.b takes the form of Eq. p2p to this accuracy. 
Similarly we obtain for the AC-current into tunnel con- 
tact a 



- 2 IwaWa I (cos z^q, , sin Va) ■ a 
An^Vaklra 



(36) 



/de , _ , 

dxa dx'^ dXcy, — u{x^^)Tig{x'^~Xc^,e)f''g{Xc^-Xa,e-uj)[f''Ga{e)-Ga{e-i^)f'']Wa{Xa,x'^^ 

m 



r 

which imphes Eq. (|lip in our hmit. 



